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We examine the validity of the weak coupling spin-susceptibility "glue" approximation (SSGA) 
in a two-dimensional Hubbard model for cuprates. For comparison, we employ the well-established 
dynamical cluster approximation (DCA) with a quantum Monte Carlo algorithm as a cluster solver. 
We compare the leading eigenvalues and corresponding eigenfunctions of the DCA and SSGA pairing 
matrices. For realistic model parameters, we find that the SSGA fails to capture the leading pairing 
symmetries seen in the DCA. Furthermore, when the SSGA is improved through the addition of a 
term with d-wave symmetry, the strength of this additional term is found to be larger than that of 
the "glue" approximation. 
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I. INTRODUCTION 

The pairing mechanism of cuprate superconductors has 
been a challenging problem since their discovery. To this 
day, among different scenarios, two stand out. First one 
is the Anderson's resonating valence bond scenario^, in 
which superconductivity is pictured as a Mott liquid of 
pairs formed by the superexchange interaction 2 . This 
is a strong-coupling approach and can predict many fea- 
tures of the cuprate phase diagram 3 -. Another is the body 
of weak-coupling approaches, including phenomenolog- 
ical models^, fluctuation exchange 5 -^ and random-phase 
approximation^, in which the pairing interaction, i.e., the 
"glue" , is mediated by the low-energy spin fluctuations. 

Recently, much effort has been devoted, both in experi- 
mental and in theoretical fronts, to find more compelling 
evidence for the spin-fluctuation-mediated pairing. On 
the experimental side, neutron scattering data show a 
prominent peak in the structure factor at the antifcrro- 
magnetic wave vector, relevant to d-wave pairing^. Using 
inelastic neutron scattering data to parametrize the ef- 
fective spin-susceptibility glue interaction, Dahm et. al£ 
find an excellent agreement between the numerically cal- 
culated features of the spectral function and the angle re- 
solved photoemission spectroscopy data. Moreover, van 
Heumen et. al. 1 ® find a correlation between the doping 
trends in the "glue spectra" , derived from optical con- 
ductivity data, and the superconducting critical temper- 
ature. 

On the theoretical side, the dynamics of this type of 
pairing has been recently investigated by numerous au- 
thors. For instance, employing an extended Hubbard 
model, Markiewicz and Bansilii argue that while mag- 
netic pairing mechanism is valid in cuprates, both high 
and low energies are relevant to pairing. Similar conclu- 
sion have been drawn using numerical calculations of the 
Hubbard model which often involve the dynamical mean- 
field treatments of the smallest system relevant to d-wave 
pairing, the cluster of four site a 12 ' 13 . However, using sim- 
ilar techniques, others argue that only the low-frequency 
part of the pairing is importan t 14 ' 15 . 

In this work, we examine the validity of the spin- 



susceptibility "glue" approximation (SSGA), expressed 
in a form similar to that of random phase approxima- 
tion, 



7,13.16-19 



tSSGA 



(K\K') = ~U 2 Xs (K' - K) 



(1) 



by exploring the momentum dependence of pairing and 
comparing it to the results obtained from a dynamical 
cluster approximation (DCA)22r— simulation. Here, V 
is the particle-particle irreducible vertex function and \s 
is the fully dressed spin-susceptibility. K = (K,cj„) de- 
notes both momentum and frequency and U is an effec- 
tive Coulomb interaction. Unlike most previous calcu- 
lations, we employ a relatively large cluster, the 16-site 
cluster, allowing for more pairing symmetries. Note that 
in our SSGA, \s is also obtained from the DCA simula- 
tion. This emulates the use of experiment to parametrize 
the "glue" approximations^. 

We find that when a finite next-nearest-neighbor 
hopping, t' , appropriate to describe the hole-doped 
cuprates 2 ^— is considered, the SSGA form of the in- 
teraction leads to p-w&ve pairing — , while for the same 
parameters the DCA yields robust d-wave pairing. We 
show that this is due to the predominant scattering with 
the antiferromagnetic wave vector, resulting from the mo- 
mentum dependence of the spin susceptibility, and the 
strong suppression of the density of states (DOS) at the 
antinodal points (strong pseudogap)22,. To re-establish 
the d-wave pairing symmetry and the agreement with 
the DCA, the SSGA can be improved by adding an addi- 
tional term with d-wave functionality in the momentum 
space. However, we find that the strength of this addi- 
tional term should be larger than that of the SSGA. 



II. FORMALISM 

We consider a two-dimensional Hubbard Hamiltonian 
H = -^2 hi ( c L c ja + h.c.) + uY^ n n n *l (2) 

ija i 

where Uj is the hopping matrix, cJ CT (cj CT ) is the creation 
(annihilation) operator for electrons on site i with spin 



2 




T/t T/t 



FIG. 1: Leading eigenvalues of the DCA [(a) and (c)] and 
SSGA [(b) and (d)] particle-particle pairing matrices versus 
temperature. Vertex function of Eq. (fjrj is used to form the 
SSGA pairing matrix. Top panels correspond to the Hubbard 
model with only nearest-neighbor hopping (t). In lower pan- 
els, a finite next-nearest-neighbor hopping is also taken into 
account. The inset of (b) shows the effective Coulomb inter- 
action used in the SSGA vertex, which is adjusted so that 
d-wave eigenvalue in SSGA is the same as its DCA counter- 
part. 

a and rij CT = c\ a Ci a . We show results for U equal to the 
bandwidth which is believed to be a realistic value for 
modeling cuprates^2r— and at filling, n — 0.95. Calcula- 
tions at different hole-dopings show that our conclusions 
are valid in the under-doped region (n > 0.85), where 
the antiferromagnetic correlations are stronger, while at 
larger dopings, the results are inconclusive since calcula- 
tions are limited by the sign problem. 

The DCA is a cluster mean-field theory that maps 
the original lattice model onto a periodic cluster of size 
N c = l? c embedded in a self-consistent host. Correlations 
up to a range L c are treated explicitly using a quantum 
Monte Carlo (QMC) solver, while those at longer length 
scales are described at the mean-field level. Previous 
DCA simulations have shown a robust d-wave supercon- 
ductivity for the Hubbard model with U comparable to 
the bandwidth. The DCA has also established pseudogap 
and antiferromagnetic phases for this model which are in 
good qualitative agreement with experimental phase di- 
agram of cuprate a 33 ' 34 . 

To study the pairing, we calculate the eigenfunctions 
of the paring matrix, Txcr^, 

£ T{K\K') Xo {K')ct>{K>) = \4>{K) (3) 

c K , 

where T is temperature, xo (= —G(K)G(—K)) is the 
particle-particle bubble in the pairing channel and V can 
be either T SSGA [Eq. [1] or calculated in the DCA. For the 
latter, we measure the two-particle Green's function in 
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FIG. 2: (Color online) The product xo(K)xo(K + Q) with 
Q = (-7T, 7r) and ui n = c*V = ttT at four different K points 
in the first Brillouin zone (1BZ) for (a) t' = and (b) t' = 
— 0.3t versus temperature. At low temperatures, the leading 
pairing symmetries in SSGA can be determined by the value 
of the product xo(K)xo(K + Q) at each K point. The one for 
K = (-7T, 0) corresponds to the d-wave symmetry and shows a 
significant decrease with t' at low temperatures. 

the pairing channel (x) in the QMC process. Then, using 
the Bethe-Salpeter equation, T is calculated by subtract- 
ing the inverse of x from the inverse of the bare bubble 
in the same channel (T — Xo 1 ~ X -1 )- The eigenfunc- 
tion <f)(K) represents the gap function and provides in- 
formation about the symmetry and the frequency depen- 
dence of the pairing. The singularities in the two-particle 
Green's function, \ = Xo/(l — Txo), which signal pair- 
ing instability, take place when the eigenvalue, A, goes to 
unity. Therefore, one can explore the superconducting 
tendencies by studying the temperature dependence of 
the leading eigenvalues. 



III. RESULTS 

Unlike DCA calculations which yield d-wave pairing, 
when a finite next nearest neighbor hopping, appropriate 
to describe hole-doped cuprates, is considered the SSGA 
pairing vertex alone [Eq. [1] does not result in pairing 
with d-wave symmetry. In Fig.[Tl we compare the leading 
eigenvalues of the SSGA and DCA pairing matrices for 
zero and finite When t' = 0, the low-temperature 
DCA results have a d-wave leading eigenvalue, followed 
by s-wave and p-wave eigenvalues. Similar results are 
found with the SSGA pairing matrix, except that s-wave 
is not one of the leading eigenvalues. However, when t' — 
— 0.3i, p-wave instead of ci-wave becomes the dominant 
pairing symmetry for the SSGA while d-wave is still the 
leading eigenvalue in the DCA. Thus, the SSGA fails to 
capture the symmetry of the pairing obtained from the 
DCA calculation. Here, U, which is shown for the two 
different values of t' in the inset of Fig[T] (b), is adjusted 
so that d-wave eigenvalue in SSGA is the same as its 
DCA counterpart. 

The momentum dependence of Xs and the renormal- 
ization of particle-particle bubble due to finite t' are re- 
sponsible for the disagreement between the symmetry of 
the leading eigenfunctions in the SSGA and the DCA. 
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FIG. 3: (Color online) The 1BZ for N c = 16. Arrows show the 
four independent (it, 7r) scatterings between K and K' on this 
cluster. The one that connects (tt, 0) to (0, n) is associated 
with the d-wave symmetry. 

To better understand why the dominant pairing in the 
SSGA is not d-wave, we employ the following approxima- 
tion: Xs{K' - K) fs Xs(Q,0)<5 K '-k,q<5(>V - u n ). This 
is motivated by the fact that the spin-susceptibility is 
considerable only at the antiferromagnetic wave vector, 
Q = K' — K = (it, tt) and at small Matsubara frequency. 
In this approximation, Eq. ([3]) can be written as 

3T 

— U 2 Xs(Q)xo(K + Q,w„)0(K + Q,u n ) « A0(K, Wn ). 

(4) 

Considering that 
3T 

— U 2 Xs (Q) X o(K,u; n )^K,Un) « A0(K + Q iUn ) (5) 

is also true, one gets 

A 2 oc Xs(Q) 2 Xo(K, w„)xo(K + Q,u) n ). (6) 

This suggests that the leading eigenvalue of the SSGA 
pairing matrix corresponds to a momentum, K, for which 
the quantity xo(K,w„)xo(K + Q,w„) has its largest 
value. Since the bubble xo(K,cj„) falls rapidly with fre- 
quency, only the lowest Matsubara frequencies are rele- 
vant for determining the leading eigenvalues. In the fol- 
lowing, we will discuss the behavior of Xo(K, w n )xo(K + 
Q,w n ) at u n = ±ttT. 

When = 0, xo(K)xo(K + Q) is the largest for 
K = (tt, 0) and K = (tt/2,tt/2) at low temperature, 
as can be seen in Fig. [5] (a). The former situation fa- 
vors d-wave pairing whereas the latter favors p-wave pair- 
ing. The close values of xo(K)xo(K + Q) at these mo- 
menta explains the competition between <i-wave and p- 
wave symmetries in SSGA. Note that for a small 2x2 
cluster, the resolution in momentum space is poor (e. g. 
K = (7r/2,7r/2) is not represented in the Brillouin zone) 
and the scattering between the nodal points which fa- 
vors p-wave pairing is suppressed. Thus, a larger cluster 
with good momentum resolution is important in captur- 
ing symmetries other than d-wave in the pairing channel. 



FIG. 4: (Color online) Fractions of the d-wave pairing inter- 
action for the SSGA term (first term) and the a term (second 
term) of Eq. for zero and finite t' versus temperature. 
Vd = V d SSGA + a is the total d-wave projected interaction 
of Eq. ([7}. For t = 0, the contribution of the a term to 
the d-wave interaction is insignificant while for t' = — Q.3t, it 
becomes more important than the SSGA term at low temper- 
atures. 



The 16-site cluster provides four independent values for 
Xo(K)xo(K + Q). This is illustrate in Fig. [3] where each 
arrow represents a Q = (tt, tt) scattering between K and 
K + Q. 

A finite t' strongly suppresses xo(K)xo(K + Q) at the 
antinodal points while it has a small effect at other k 
points [see Fig. [H(b)]. Therefore, according to Eq. ©, 
the d-wave eigenvalue will also be suppressed relative to 
the p-wave eigenvalue. Consequently, the p-wave pairing 
will be dominant in the SSGA, which explains the results 
discussed in Fig [1] The renormalization of the bare bub- 
ble at the antinodal points can be understood from the 
changes in the DOS. At low temperatures, the low energy 
DOS at K = (0, 7r) is strongly suppressed [i.e., the pseu- 
dogap is enhanced] with t' for the hole-doped systems 
while this effect is negligible at other k points 29 . This in- 
dicates that in the over-doped region, where pseudogap is 
less pronounced, SSGA might be a good approximation. 

We find that the spin-susceptibility representation of 
the pairing interaction, which is appreciable only at 
K' K = (7r,7r), does not always yield d-wave as the 
leading pairing symmetry. This seems to be especially 
true when realistic parameters for cuprates such as next- 
nearest-neighbor hopping are considered. We show that 
despite the inherent correspondence between the Q = 
(tt, tt) wave vector, at which the SSGA interaction is 
large, and the d-wave symmetry, other pairing symme- 
tries which also associate with the (tt, tt) wave vector can 
be dominant. On the other hand, the DCA exhibits a 
robust d-wave pairing for the same physical parameters 
of the Hubbard model. 

To investigate the possibility of a missing term, we pro- 
pose to add an extra term to r SSGA which enhances the 
scattering between the antinodal points. Such approxi- 
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mation for the interaction vertex can be written as: 

T'(K\K') = ^U 2 Xs (K' -K)- a<t> d {K')<j) d (K). (7) 

Here, <f>d is the d-wave eigenfunction of the DCA pairing 
matrix, and U and a are temperature-dependent fitting 
parameters which are adjusted to reproduce both the d- 
wave and the p-wave leading eigenvalues of the DCA pair- 
ing matrix. It is worth mentioning that Maier et. al*£ 
proposed a similar term, — Jg(K)g(K'), to be added to 
the SSGA form in which g(K) oc (cosK x — cosK y ) is the 
d-wave form factor and J is an effective exchange interac- 
tion. However, their motivation for the necessity of this 
extra term is quite different from ours. Since they did 
not consider a finite t' , the SSGA gave a d-wave pairing. 
So, this term was suggested only to restore the large fre- 
quency behavior of the pairing gap found in the DCA. 
Our focus, on the other hand, is to examine the rele- 
vance of an additional term based on a momentum space 
argument. Note that because our calculations are done 
using QMC on the imaginary frequency axis, generally 
comments about the frequency dependence of the pair- 
ing interaction cannot be made. However, if relevant, the 
additional term may be responsible for the instantaneous 
part of the interaction, as described in Ref. [IH We also 
find that using (f>d(K) instead of g(K) in our approxi- 
mation provides a better fit to the DCA results while 
imposing the <i-wave symmetry. 

For a finite t 1 , the second term in Eq. (|7|) (a term) has 
a prominent role in capturing the correct leading sym- 
metries in this new form of the interaction. In Fig. 0] we 
show the fractional values of a and the d-wave component 
of the SSGA texmi&36, 

ySSGA = _ £ M K<)^U\ S (K> - K)MK), (8) 

K,K< 

for the two values of t' . We define Vd = V c f SGA + a as 
the total d-wave projected interaction. When t' = 0, the 



contribution of the a term to the d-wave interaction is 
insignificant. However, when t' — — 0.3i, this term has 
a dominant role in the d-wave pairing at low tempera- 
tures. So, although the additional term enhances d-wave 
pairing, its contribution overshadows the contribution of 
the main part of the interaction i.e. the term propor- 
tional to the spin-susceptibility. This suggests that the 
SSGA may not be improved by simply adding terms that 
enhance d-wave scattering. 

IV. CONCLUSIONS 

We study the validity of the SSGA representation of 
the pairing interaction in the Hubbard model. By com- 
paring the leading pairing symmetries of this interac- 
tion with the pairing symmetries produced by the DCA, 
we find that this approximation alone does not capture 
the correct pairing symmetry in the under-doped and 
optimally-doped regions, particularly when a finite t' is 
considered. We do not dismiss the possibility that SSGA 
can be valid in the over-doped region, where the pseu- 
dogap phenomenon is less relevant. We show that this 
form of the interaction, which is large only at the antifcr- 
romagnetic wave vector, requires an additional term with 
d-wave symmetry to yield d-wave pairing at low temper- 
atures. However, in case of finite t' , the additional term 
dominates the interaction as the temperature is lowered. 
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